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Light-mediated silencing and stimulation of cardiac excitability, an important complement to electrical 
stimulation, promises important discoveries and therapies. To date, cardiac optogenetics has been studied 
with patch-damp, multielectrode arrays, video microscopy, and an all-optical system measuring calcium 
transients. The future lies in achieving simultaneous optical acquisition of excitability signals and 
optogenetic control, both with high spatio-temporal resolution. Here, we make progress by combining 
optical mapping of action potentials with concurrent activation of channelrhodopsin-2 (ChR2) or 
halorhodopsin (eNpHR3.0), via an all-optical system applied to monolayers of neonatal rat ventricular 
myocytes (NRVM). Additionally, we explore the capability of ChR2 and eNpHR3.0 to shape 
action-potential waveforms, potentially aiding the study of short/long QT syndromes that result from 
abnormal changes in action potential duration (APD). These results show the promise of an all-optical 
system to acquire action potentials with precise temporal optogenetics control, achieving a long-sought 
flexibility beyond the means of conventional electrical stimulation. 

Optogenetics has revolutionized neuroscience' '', offering exquisite control of cell excitability with a 
resolution and flexibility complementary to more traditional strategies involving electrical stimulation. 
In particular, optogenetics affords precise spatiotemporal control via light-activated unipolar injection 
of depolarizing''"" (ChR2) and hyperpolarizing currents^''' (eNpHR3.0), kindling strong interest in applications 
beyond neurons, such as in the heart. Hoffmann et al.' demonstrated light-induced stimulation of confluent 
layers of HL-1 cells expressing ChR2, recording electrical propagation with a microelectrode array. Bruegmann et 
al.*' undertook light-induced stimulation of ChR2-expressing murine muscle and heart. Arrenberg et al.'^ trans- 
genically expressed NpHR and ChR2 H134R in zebrafish cardiomyocytes to locate the pacemaker in the devel- 
oping heart, and to simulate pathological cardiac conditions like tachycardia, bradycardia, atrioventricular block, 
and cardiac arrest. AbUez et al. recorded photostimulation and photoinhibition of cardiomyocytes derived from 
human embryonic stem cells'" and human induced pluripotent stem cells". Jia et al.'^ used an all-optics system to 
acquire optical recordings of Ca-* transients from tandem units of cardiomyocytes with HEK293 cells expressing 
ChR2. Nussinovitch et al.''' co-cultured cardiomyocytes with NIH-3T3 fibroblasts expressing ChR2, and used 
microelectrode arrays to monitor resynchronized electrical activity after conduction block. Kirkton et al.'* 
activated engineered HEK293 cells expressing ChR2 in a cardiac setting. Finally, much interest lies in the 
computational modeling of optogenetics control'" '"'". 

In all these endeavours, an aU-optics system would be ideal, combining the strength of optogenetics (precise 
spatiotemporal control and unipolar current injection) and optical mapping of excitability signals'" (noninvasive, 
high-resolution visualization of conduction patterns). Thus far, however, only Ca"^* transients have been mea- 
sured in an all-optics manner'^", whereas optical monitoring of the more fundamental voltage signals that 
comprise action potentials has yet to be reported for this configuration. Here, we therefore devise an apparatus 
that permits readouts of action-potential morphology and conduction, concurrent with orthogonal activation of 
optogenetic constructs expressed in NRVM monolayers. As an initial application, we explore the capability of this 
all-optics system to shape cardiac action-potential morphology, in particular to lengthen and shorten APDs. Such 
a capability might ultimately be desired to probe and manipulate pathologic conditions relating to APD abnor- 
malities. In particular, long and short QT syndromes (LQTS or SQTS) derive respectively from unusually 
prolonged or abbreviated APDs, predisposing for conduction heterogeneities, ventricular arrhythmias, and 
sudden cardiac death^""^^. The efficacy of current treatments varies considerably with underlying pathogenesis 
and severity of arrhythmias^', and further insight into the underlying mechanisms, spanning spatiotemporal 
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scales from molecules to cell networks, would be desirable. Our all- 
optics system may furnish helpful tools towards this long-term goal. 

Results 

Biophysical properties of ChR2 and eNpHR3.0 during action- 
potential waveforms. To assess the overall feasibility of using opto- 
genetic constructs to influence cardiac action-potential morphology, 
we initially assessed the magnitude of optogenetic currents using 
whole-cell patch-clamp of HEK293 cells stably expressing either 
ChR2 or eNpHR3.0. For ChR2 (Fig. la), we transiently activated 
current with pulsed blue-light illumination (447 nm) while voltage 
was held constant. In different trials, voltage was set at different levels 
ranging from —80 mV to 80 mV in 10 mV steps. Activation of 
inward current was observed between —80 mV to 10 mV, and 
current amplitude decreased as voltage became more positive. 
Outward current was present from 20 mV to 80 mV, and current 
amplitude increased as voltage became more positive. These data are 
consistent with a reversal potential between 10-20 mV. In agree- 
ment with previously reported patch-clamp results^, partial inactiva- 
tion of currents was clearly present. For eNpHRS.O (Fig. lb), currents 
were transiently activated by pulsed green-light illumination 
(530 nm) delivered within an analogous protocol. Light activation 
produced relatively uniform inward chloride currents at all voltages. 
More detailed inspection reveals that the chloride current increased 
slightly with increasing voltage, indicating that eNpHR3.0 is weakly 
voltage-dependent as expected from its pump-like properties. In 
agreement with previously reported patch-clamp results'", there 
was no inactivation of current. Population data for ChR2 currents 
are shown for peak (Fig. Ic, red) and steady-state (Fig. Ic, black) 
current levels; these illustrate the inwardly rectifying character of this 
channel under physiological ionic gradients''"''. As well, population 
data for eNpHR3.0 currents are shown for peak (Fig. Id, red) and 
steady-state (Fig. Id, black) current levels"'; these indicate the near 
step-like kinetics of this pump. The biophysical properties were also 
characterized under variable light intensities to show that net current 
increases with increasing light intensity for both ChR2 and 
eNpHR3.0 (Supplementary Fig. S2). Additionally, the inactivation 
kinetics of ChR2 show increasing inactivation with increasing 
voltage or blue-light intensity (Supplementary Fig. S3). 

To further gauge the potential suitability of the optogenetic con- 
structs for altering action- potential morphology, currents were eli- 
cited during voltage-clamp with cardiac action-potential waveforms 
(APW), as shown in Figures le and If Optogenetic currents were 
allowed to reach steady state with respect to light exposure prior to 
trace onset (Figs, le-f, note 500-ms axis breaks in light-pulse wave- 
form). As expected, based on the voltage step recordings, ChR2 
current was greater at negative voltages (Fig. le) and current 
decreased with increasingly depolarized voltages. The overall current 
trace morphology was similar to the APW. Similar to the constant 
voltage step recordings (Fig. lb), eNpHR3.0 current acquired during 
an APW (Fig. If) was also relatively constant, indicating little volt- 
age-dependence of current. 

Optical mapping of action potentials with orthogonal ChR2 or 
eNpHR3.0 excitation. A significant challenge was to record trans- 
membrane potential via a voltage-sensitive fluorescent dye, while 
independently switching on/off ChR2 or eNpHR3.0. Specifically, 
the absorption spectrum of the more commonly used voltage- 
sensitive dye, di-4-ANEPPS^'', fully overlaps the absorption spec- 
trum of ChR2', and partially overlaps the absorption spectrum of 
eNpHR3.0^. This overlap makes it impossible to independently excite 
di-4-ANEPPS and the optogenetic constructs. The solution was to 
use the long-wavelength voltage-sensitive dye PGHP"" instead, with 
an absorption spectrum"" separate from that of ChR2 and eNpHR3.0 
(Supplementary Fig. S4). 
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Figure 1 | Biophysical properties of optogenetic constructs and their 
suitability for altering excitability during action potentials. ChR2 (a, 
c) and eNpHRS.O (b, d) were initially characterized at various 
constant levels of voltage, using whole-cell patch clamp in HEK293 
cells, (a) ChR2 supports light-gated, voltage-dependent currents with 
a reversal potential of —17 mV according to the curve fit of the 1-V 
relation in (c). Partial inactivation of currents is clearly present, 
(b) eNpHR3.0 is light-gated, but only weakly voltage-dependent, (c) 
ChR2 population data for peak (red) and steady-state (black) current 
levels, n refers to number of cells; bars show sem. (d) eNpHR3.0 
population data for peak (red) and steady-state (black) currents, n 
refers to number of cells; bars show sem. Panels (e) and (f) examine 
optogenetic currents during voltage-clamp with cardiac action- 
potential waveforms gauged from neonatal rat ventricular myocytes 
probed with PGHl as the voltage-sensitive dye. (e) ChR2 currents 
during APWs are large and enhanced by repolarization. Current 
amplitude between the maximum during the trace onset at —80 mV, 
and the minimum near the peak of the APW, was 427 ±159 pA 
(mean ± sem, n = 5 cells), (f) eNpHR3.0 currents during APWs are 
smaller but largely independent of voltage. Current amplitude at the 
trace onset 32 ± 0.4 pA (mean ± sem, n = 12 cells). All blue 
illumination at 7.2 mW/mm". All green illumination at 1.4 mW/mm^. 
Axis breaks (500 ms) in light-power schematics atop (e) and (f), 
ensure that channels had reached near steady-state (in)activation with 
respect to light, prior to trace onset. 
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Figure 2 | Optical mapping system featuring detection of red-shifted 
PGHl voltage-sensitive dye signals, concurrent with orthogonal 
excitation of either ChR2 or eNpHR3.0 optogenetic constructs. 



(a) Optical layout with dual-pathway excitation modules placed above 
NRVM monolayer. Fluorescence emission is captured from below, 
through a 253-channel fiber optic bundle, (b) Top, Exemplar NRVM 
(without optogenetic expression) action potential averaged from multiple 
channels, using customary di-4-ANEPPS voltage-sensitive dye. Bottom, 
corresponding activation map of uniform propagation. The color bar 
represents activation time in ms defined as the time of the stimulus to 50% 
of peak, (c) Top, Action-potential morphology was largely maintained 
with use of the red-shifted PGHl dye. APDs were modestly prolonged 
compared to di-4-ANEPPS. Bottom, Propagation with PGHl was also 
uniform throughout the monolayer, though PGHl -stained monolayers 
exhibited somewhat slower conduction, (d, e) Blue- or green-light pulses 
had no collateral effects on NRVMs lacking optogenetic expression, as 
demonstrated by exemplar PGHl traces from a single recording channel. 
Blue or green illumination caused baseline shifts of the fluorescence signal, 
requiring crosstalk compensation, which is described in more detail in 
Method "Gap compensation for blue- and green-light bleed-through." 
AH blue illumination at 7.2 mW/mm^. All green illumination at 
1.4 mW/mm^ 

An optical mapping system'" was engineered (Fig. 2a) to enable 
detection of red-shifted PGHl dye signals"^^^ concurrent with ortho- 
gonal excitation of either ChR2 or eNpHR3.0. NRVM monolayers 
stained with PGHl were directly illuminated using a red (655 nm) 
LED array passed through an interference filter (690/60 nm). Action 
potentials, in the form of fluorescence emission, were acquired at 
wavelengths above 760 nm using a long- pass filter and a 253-channel 
fiber-optic bundle. To excite ChR2, light from a perpendicularly 
aligned blue (447 nm) LED array was passed through an interference 
filter (440/60 nm), and directed onto the monolayer using a dichroic 
mirror splitting at 475 nm. To excite eNpHRS.O, light from a per- 
pendicularly aligned green (530 nm) LED array was passed through 
an interference filter (560/50 nm), and directed onto the monolayer 
using a dichroic mirror splitting at 618 nm. The resulting data fea- 
tures multiple channels of action -potential data. Each channel 
reports the signal of a single fiber-optic, yielding a composite action 
potential from many syncytial myocytes in a 1-mm diameter spot of 
the monolayer. 

To verify that PGHl does not substantially affect action-potential 
morphology, spatially averaged action potentials acquired with the 
traditionally used di-4- ANEPPS (Figs. 2b top and 2c top, black) were 
compared to those obtained with PGHl (Fig. 2c top, red). Overall, 
action potentials measured with PGHl were similar to those 
observed with di-4-ANEPPS, with PGHl responses exhibiting 
slightly longer APD and modestly slower depolarization/repolariza- 
tion rates. Activation isochrone maps with di-4-ANEPPS (Fig. 2b 
bottom) and PGHl (Fig. 2c bottom) both indicated uniform pro- 
pagation, with possibly slower conduction velocity measured using 
PGHl. This slight slowing of action-potential upstroke and down- 
stroke observed with PGHl, along with a somewhat diminished 
conduction velocity, may relate to the higher levels of PGHl staining 
required to achieve an adequate signal-to-noise ratio. Nonetheless, 
these effects were small enough to proceed with the use of PGHl in 
the experiments that foUow. 

To confirm that blue or green illumination for optogenetic stimu- 
lation would not in itself interfere with action-potential responses 
recorded by PGHl, NRVM monolayers lacking optogenetic con- 
structs were stained with PGHl and electrically paced. Action poten- 
tials were then recorded without (Fig. 2d left) and with (right) 
delivery of multiple blue- or green-light pulses. As baseline, the con- 
trol response on the left (black trace) was fitted with the smooth red 
function for convenient comparison. Multiple high-intensity blue- 
light pulses introduced bleed-through artifacts manifest as gaps in 
the black trace on the right (see Methods "Gap compensation 
and Supplementary Fig. SI). Nonetheless, superposition of the 
smooth red curve fit from control indicates that the blue pulses did 
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Figure 3 | Transduction of NRVM monolayers with lentiviruses bearing 
ChR2 or eNpHRS.O did not appreciably perturb action-potential 
morphology from control (absent optogenetic stimulation with blue or 
green light), (a, d) Control NRVMs without viral transduction, 
(a) Transmission (top) and fluorescence (bottom) confocal images. 

(d) Action potentials were averaged from multiple channels of optical 
mapping data, with PGHl as dye. Parallel culture to those transduced with 
lentiviruses in later panels, (b, e) Results for NRVM transduced with 
lentivirus bearing ChR2; format as in (a), (b) Fluorescence is from YFP 
fused to ChR2, and reflects expected membrane targeting of ChR2. 

(e) Action potential (blue) essentially unchanged from control (gray), 
(c, f) Results for NRVM transduced with lentivirus bearing eNpHR3.0; 
format as in (a), (c) Fluorescence is from YFP fused to eNpHR3.0, and 
reflects membrane targeting of optogenetic construct, (f) Action potential 
(green) essentially unchanged from control (gray), as copied from panel 
(a). 

not by themselves affect action-potential morphology. Likewise, 
multiple high-intensity green-light pulses did not by themselves 
affect action-potential morphology (Fig. 2e). 

Un-stimulated ChR2 or eNpHR3.0 minimally alters cardiac action 
potentials. To exclude artifactual changes in action-potential mor- 
phology by expression of optogenetic constructs in the absence of 
their activation, we compared NRVM monolayers in the absence of 
lentiviral transduction (Fig. 3a) with those transduced with lentivirus 
bearing ChR2 (Fig. 3b) or eNpHR3.0 (Fig. 3c). Under optical map- 
ping using PGHl, averaged action potentials for the un-transduced 
monolayers displayed their usual morphology as shown in Figure 3d. 
Throughout, no blue or green illumination was applied to exclude 
activation of optogenetic constructs. Infection of NRVM monolayers 
with lentivirus bearing ChR2 resulted in widespread expression 
across cells, each showing the expected membrane targeting of YFP- 
tagged channels (Fig. 3b, bottom panel). The spatially averaged 
action potential (Fig. 3e, blue) from the same culture and infection 
indicates no appreciable perturbation of action-potential morpho- 
logy from control (Fig. 3e, gray trace). Similarly, infection of NRVM 
monolayers with lentivirus bearing eNpHR3.0 was efficient and 
showed membrane targeting of the YFP-tagged pumps (Fig. 3c, 
bottom panel). The spatially averaged action potential (Fig. 3f, 



green) also suggests minimal perturbation of action-potential mor- 
phology from control (Fig. 3f, gray trace, reproduced from Fig. 3d). 
Infection with empty virus (only GFP) also did not affect action- 
potential morphology (Supplementary Fig. S5). Thus, expression of 
optogenetic constructs did not in themselves alter action-potential 
behavior. 

Functional excitation or silencing of NRVM monolayers with 
ChR2 or eNpHR3.0. We initially tested for simple stimulation and 
silencing of action potentials by brief activation of the optogenetic 
constructs. Voltage was monitored with optical mapping utilizing 
PGHl. In monolayers transduced with ChR2, action potentials 
were elicited by a train of blue-light pulses (51 ms duration) at 
maximal power (7.2 mW/mm^) with no electrical stimulation 
(Fig. 4a, left). The intensity and duration of the blue-light pulses 
were long enough to simultaneously excite the entire monolayer, 
preventing propagation (Fig. 4a, right). Accompanying manipula- 
tion of monolayer contraction was also observed (Supplementary 
Movie 1 shows that a monolayer bearing ChR2 contracts after 
short-duration, high-intensity blue-light pulses). When the blue- 
light pulse intensity was lowered or the duration was shortened, 
propagation was evident across the monolayer (data not shown). 
In eNpHR3.0-expressing NRVMs, electrically paced action 
potentials were silenced with variable intensity green-light pulses 
(Fig. 4b, top). The green-light intensities (Fig. 4b, green trace, in 
mW/mm") were 1.4, 1.4, 1.24, 0.67, 0.14, 1.4, and 1.4. Electrical 
stimulation nicely induced a propagating depolarizing wave from 
left to right (Fig. 4b, bottom, leftmost and rightmost optical maps). 
The first four green-light pulses successfully abolished the electrically 
paced action potentials, also evident across space by the lack of 
activation in the isochrone map (Fig. 4b, bottom, optical map 
second from left). However, with each of these four successful 
sUencing events, an action potential was elicited shortly after 
cessation of the green-light pulse, much like an anode-break-like 
response. As might be expected, this break-like response produced 
nearly uniform excitation across the preparation (Fig. 4b, bottom, 
map third from the left). Analogous manipulation of monolayer 
contraction was also observed (Supplementary Movie 2). The fifth, 
low intensity green-light pulse failed to silence the electrically paced 
action potential, hinting at the possibility of tuning green-light 
intensity to do more than simply switch off electrically stimulated 
responses, an issue that will be explored two sections below. 

APD elongation using multiple blue-light pulses in ChR2-expressing 
monolayers. Up to this point, we have demonstrated the ability to 
record action potentials under an optical mapping configuration 
permitting orthogonal excitation of optogenetic constructs, repro- 
ducing previously reported optogenetic phenomena in cardiac 
preparations. Here, we directly test whether ChR2 may be used to 
sculpt action-potential morphology itself We sought to increase 
APD using blue-light excitation in ChR2-expressing monolayers. 
First, to verify that APDs can be increased appreciably in NRVMs, 
and to set a heuristic target for optogenetic APD elongation, we 
performed optical mapping with PGHl on monolayers incubated 
in 1 |.iM Bay K 8644, an L-type calcium channel agonist"^'. Figure 5a 
indicates substantial augmentation of APD, as illustrated by com- 
parison of the spatially averaged action potential obtained under 
pharmacological enhancement of L-type current (Fig. 5a, right, 
black trace) to that of control action potentials (Fig. 5a left, 
and right gray trace). Conversely, a calcium channel antagonist 
shortened APD (Supplementary Fig. S6). To test for elongation of 
APD by ChR2, increasing numbers of punctuated blue-light pulses 
(7.2 mW/mm^) were administered during the repolarization phase 
of electrically paced action potentials in a monolayer expressing 
ChR2 (Fig. 5b, i-iv). Responses without blue-light stimulation (o) 
were interspersed to confirm the return to a stable baseline between 
presumed perturbations. Gaps indicate periods of blue-light bleed- 
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Figure 4 | ChR2 and eNpHRS.O can respectively trigger and silence 
action potentials, as seen from optical mapping data obtained with 
PGHl. (a) Exemplar result of NRVM monolayer expressing ChR2 exposed 
to high-intensity blue pulses (51 ms, top, blue trace), without electrical 
stimulation. Hardly noticeable gaps in records reflect periods of blue-light 
bleed- through. Excitation of the monolayer occurs simultaneously with no 
propagation. Time zero in the activation isochrone map corresponds to the 
start of the blue pulse. Position of the exemplar trace is marked with a white 
asterisk on the isochrone map. (b) In eNpHR3.0-expressing NRVMs, 
electrically paced action potentials can be silenced with variable intensity 
green pulses (top, green trace). The first four green high-intensity pulses 
silence action potentials but evoke anode-break-like action potentials 
upon cessation of green pulses. Electrical pacing in gray (just below 
fluorescence trace). Bottom, isochrone maps corresponding to various 
phases of trace as shown. Propagated depolarizing wave for electrically 
stimulated responses (leftmost and rightmost maps). No excitation across 
monolayer during the second silencing green pulse (second from left map) . 
Non-propagating uniform excitation in break-like response (third from 
left map). Zero time of maps is instant of electrical stimulus, except for 
break-like response where cessation of green-light pulse is taken as time 
zero. Position of recording site for exemplar trace is marked with a white 
asterisk on the isochrone maps. 

through (see Methods "Gap compensation ..."). Even on the 
comparatively compressed timescale of Fig. 5b, APD prolongation 
during blue-light stimulation is evident. To confirm this impression. 
Fig. 5c aligns the various APWs on an expanded time-base, and fits 
to these waveforms (smooth red curves) are overlaid to improve 
visual comparison. To improve precision, fits were performed 
sequentially on trials with increasing number of blue pulses (see 
Methods "Gap compensation ..."). Seen in this manner, each 
additional blue-light pulse induced further augmentation of APD, 
as shown by the shaded red region denoting the difference between 
blue-light-stimulated waveforms and control. Alignment of smooth 
fits to the various responses, shown at the bottom right, underscores 
the graded manner in which high-intensity blue-light pulses can 



elongate APD in a graded fashion. Indeed, the waveform observed 
with four blue-light pulses (Fig. 5d, red trace) approximates the 
elongation produced by our heuristic maneuver with 1 [iM Bay K 
8644 (black trace). These exemplar trends were fuUy corroborated by 
APDgo population data from multiple optical channels on a single 
coverslip (Fig. 5e), and by APDgo data averaged over multiple 
coverslips (Fig. 5f). Notably, the APDgo measurements were made 
on directly present portions of the data that did not require smooth- 
function fits. Thus, these latter analyses (Figs. 5e-f) thereby confirm 
the ability to optogeneticaUy sculpt APD elongation by delivery of 
blue pulses. 

APD shortening using green-light pulses in eNpHR3.0-expressing 
monolayers. To explore the counterpoint of the previous experi- 
ments, this section tests for the possible shortening of action 
potentials by activation of eNpHRS.O in NRVM monolayers. 
In this regard, while high-intensity green light could activate 
eNpHR3.0 to such an extent that action potentials were altogether 
silenced (Fig. 4b), we wondered whether lower-intensity green-light 
activation might allow for responses with attenuated duration. 
Accordingly, we scrutinized electrically stimulated responses of 
NRVMs, particularly during modest activation of eNpHR3.0 
(Fig. 6a). Green-light intensities in this panel were 1.4, 1.4, 1.24, 
0.67, 0.14, and 1.4 (mW/mm^). Similar to Fig. 4b, higher intensity 
green-light pulses altogether abolished action potentials. However, 
during the weaker green-light pulses (0.67 and 0.14 mW/mm^), 
miniature action potentials with blunted amplitude and shortened 
durations were observed (Fig. 6a, green asterisk). Figure 6b displays a 
close-up comparison of a control action potential (Fig. 6a, black 
asterisk) and a miniature action potential (Fig. 6a, green asterisk). 
Clearly evident are the action potential morphology changes 
wrought by activating eNpHR3.0 with green light: a diminished 
amplitude and shortened duration at the 80% level (red dashed 
line). Population data for APDs averaged across multiple channels 
on a single coverslip (Fig. 6c) verified that the shortened and blunted 
action-potential morphology, as induced by modest green light, was 
present across the monolayer. Curiously, if green pulses were 
selectively delivered during the repolarizing phase of action 
potentials, as was done with ChR2 in Figure 5b, action potentials 
could not be shortened. Presumably, the expression of eNpHR3.0 
current was insufficient to perturb action-potential trajectories once 
native voltage-activated channels become fully activated by larger 
depolarizations. Nonetheless, these results indicate the potential 
for optogenetic hyperpolarizing currents to induce shorter action 
potentials, particularly if larger eNpHR3.0 currents could be 
produced. 

Discussion 

Our work achieves a significant technical advance in cardiac opto- 
genetics-the ability to optically map voltage signals underlying 
action-potential and conduction profiles in cardiac monolayers, 
while simultaneously modulating the observed action potentials via 
optogenetics. In particular, using an aU-optics setup, we show that 
cardiac APD can be increased in a graded manner by activating ChR2 
currents, and action-potential amplitude and duration can be down- 
wardly modulated by activating eNpHR3.0 currents. By contrast, the 
bulk of prior optogenetic effort in the heart has focused on outright 
excitation or silencing of activity^ '" '"^. As well, when measured, con- 
duction profiles were obtained with a micro-electrode array' or by 
optical mapping of calcium transients". Instead, our work here 
demonstrates that it is possible to acquire transmembrane potentials 
in this context, promising future opportunities to combine simultan- 
eous recordings of voltage and Ca"'* with orthogonal excitation of 
ChR2 and eNpHR3.0. 

One potential limitation might concern the use of ChR2 to pro- 
long APDs in preparations where the plateau phase (phase 2) is more 
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Figure 5 | APD can be increased in ChR2-transduced NRVMs, using multiple blue pulses (20 ms on, 30 ms off) delivered during the repolarization 
phase, (a) Action potential of control NRVMs (left black trace, without ChR2) can be elongated by 1 |j.M Bay K 8644 (right black trace; control 
reproduced as gray trace), demonstrating pharmacological APD prolongation, (b) Exemplar trace demonstrating APD lengthening by transient ChR2 
activation in NRVMs. Action potentials with increasing numbers of blue pulses (i-iv, blue bars) interspersed with control APs (o). Trace shows gaps due 
to crosstalk compensation of the baseline shifts during blue illumination, (c) Expanded time-base display of records in (b). APD increases in a graded 
manner with number of blue pulses. AH action potentials (control (o) and blue pulse action potentials (i-iv)) were fitted (red) to facilitate visualization 
after crosstalk compensation of blue-light pulses. Right bottom, comparison of fits, (d) APD with four blue pulses (iv) is significantly increased to 
resemble 1 |iM Bay K 8644 action potential, (e) Bar graphs of average APDgo for various types of action potentials, confirming a graded increase in APD 
with increasing blue-pulse number (mean ± sem, n = 16 channels from single coverslip). The four-blue-pulse APDgo approximates pharmacologically 
lengthened APDgo observed with 1 [iM Bay K 8644 (mean ± sem, n = 82 channels). Analogous results for Bay K 8644 were obtained in a total of /j = 9 
coverslips. Statistical significance computed via 1-way ANOVA with Tukey-Kramer multicomparison criterion, with ** signifying p < 0.05 and 
* signifying p < 0.01. (f) Similar APD trends were confimed for population analysis over n = 5 coverslips (parentheses), where APDgo values were 
normalized to values obtained from responses without a blue-light pulse. Statistical significance computed via 1-way ANOVA with Tukey-Kramer 
multicomparison criterion, with ** signifying p < 0.05 and * signifying p < 0.01. 



SCIENTIFIC REPORTS | 4:6125 | DOI: 1 0. 1 038/srep061 25 



6 



a Green pulses 




pulse pulse like pulse 
APs APs APs APs 

o o o 

Figure 6 | Electrically stimulated action potentials become shorter and 
smaller in amplitude during low-intensity green pulses delivered to 
NRVMs expressing eNpHR3.0. (a) Effects of variable intensity green 
pulses (top). Electrical stimuli (bottom). Reference 'pre-pulse' response 
(black asterisk) is electrically paced. Higher intensity green pulses (first 
two) silence action potentials. Lower intensity pulse (green asterisk) yields 
a short duration and smaller amplitude action potential. Anode-break-like 
response (red asterisk). Post-green pulse (blue asterisk) shows return to 
control, (b) Closer view of control (black asterisk) compared to low 
intensity green-pulse action potential (green asterisk), (c) APD metrics 
(mean ± sem, n = 12 channels from single coverslip), confirming APD 
shortening by low intensity green-light pulses (green asterisk). Green pulse 
APD (green asterisk) calculated at APDgo of control ((b), red dashed line) 
without normalization. Statistical significance computed via 1 -way 
ANOVA vnth Tukey-Kramer multicomparison criterion, with * signifying 
p < 0.01. n.s. signifies p > 0.25. 

pronounced (e.g., humans, rabbits, and guinea pigs) than in the rat 
NRVMs used here (similar to mouse). Because ChR2 current 
reverses at about 10-20 mV (Figure IC), one might wonder whether 
the resulting small magnitude of ChR2 currents during a human-like 
plateau phase (near 10-20 mV) would preclude appreciable APD 
prolongation by ChR2 in this setting. However, it is worth consider- 
ing how 7ki> which manifests diminutive currents near its reversal 
potential at £k. nonetheless anchors the cellular resting membrane 
potential near its reversal potential. Moreover, the magnitude of 
native currents and conducatance during the plateau phase of 
human-like action potentials would be relatively smalP", so ChR2 
need overcome only moderately sized native currents to prolong the 
plateau phase of these action potentials. Favoring this outcome on 
empirical grounds is the ability of the 'first blue light pulse' in 
Figure 5C to prolong the NRVM action potential. In particular, the 
absolute potential at this point of the action potential would be 
greater than 0 mV^''™. Thus, we would anticipate that ChR2 could 



still prove effective at manipulating action potentials in preparations 
manifesting a more pronounced plateau phase. 

A more stringent limitation concerns practical limits on the levels 
of eNpHR3.0 current that could be expressed in NRVMs. In particu- 
lar, action potential duration and amplitude could be attenuated only 
by commencing modest activation of eNpHR3.0 prior to an electrical 
stimulus pulse, and continuing this activation throughout a sub- 
sequent action potential. Notably, we were unsuccessful in altering 
action potential shape by introducing green-light pulses only during 
the plateau or repolarizing phases of action potential, as was accomp- 
lished with ChR2. This curious dichotomy likely arises because the 
magnitude of hyperpolarizing photocurrent required to blunt action- 
potential waveforms is less if action potentials are attenuated prior to 
their fuU evolution. Delivery of modest hyperpolarizing photo- 
current prior to an electrical stimulus may allow the voltage thresh- 
old for eliciting an action potential to be reached, but this delivery 
may stUl limit the ensuing extent of voltage activation of /Na and /ca,L 
currents, which in turn influences action-potential morphology. By 
contrast, once a 'normal' action potential reaches peak amplitude, 
comparatively larger and 7ca,L currents must be dealt with, and 
somewhat larger repolarizing /ks and /kt currents may dwarf the 
relative contribution of repolarizing photocurrents achievable by 
eNpHR3.0 (Appendix B). Unfortunately, we could not produce lar- 
ger optogenetic hyperpolarizing currents, because eNpHRB.O cur- 
rents were smaller to start (cf. Fig. If in HEK293 cells), and 
eNpHR3.0 also seemed more toxic to myocytes than was ChR2. 
An alternate hyperpolarizing photocurrent, Archaerhodopsin", 
was tested but currents were yet smaller in NRVMs. 

Another current challenge regards the finite spectral range within 
which the dyes and optogenetic constructs must function. For the 
near future, the issues of spectral crosstalk wiU be resolved by the 
development of better performing voltage-sensitive dyes with 
absorption spectra at wavelengths outside the blue and green 
range^''"'''. Improved genetically-encoded voltage indicators are also 
now becoming available''' Furthermore, the toolbox of optoge- 
netic constructs continues to grow, diversifying the types of light- 
activated control available for basic science research"'. As for the 
limitations of eNpHR3.0 current amplitudes, one simple improve- 
ment would be to deliver greater green-light activation; the max- 
imum green-light power we were able to achieve via commercial 
LED components was —1.4 mW/mm^, and our characterization of 
eNpHR3.0 indicates that further eNpHR3.0 current could be 
recruited with greater green-light excitation (Supplementary Fig. 
S2). Another potential improvement would be to use transgenic mice 
or rats bearing the optogenetic constructs, thus offering poten- 
tially uniform and highly expressing optogenetic constructs across 
myocytes. 

More broadly, multiple exciting applications loom on the horizon. 
Firstly, with regard to arrhythmia research, a next-generation version 
of our apparatus using digital micromirror devices' could create and 
study customized time-varying spatial patterns to investigate how 
local geometry interacts with cellular excitability to promote or 
inhibit reentrant waves. The flexibility and power of such a strategy 
are enormous. Secondly, another relevant dimension of arrhythmo- 
genesis concerns how many cells must manifest ectopy or altered QT 
properties, and in what distribution, to induce appreciable arrhyth- 
mogenic behavior at the whole-heart and ceU-monolayer level. This 
question has received much attention from the computational bio- 
logy community" ""'^ but means of experimental verification are at 
present limited"'*'. Optogenetic approaches to trigger ectopy or 
altered action-potential morphology in variable numbers of cells, 
with differing proximity and distribution, are likely to weigh impor- 
tantly, much as optogenetics has commented on the minimal num- 
ber of neurons required to form memories™. Thirdly, ChR2 promises 
precision prolongation at certain points of the action potential wave- 
form, a feature difficult to achieve by injection of external currents, or 
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by manipulating the complement or behavior of ion channels. This 
capability is particularly well suited to testing for the existence of 
critical windows for inducing proarrhythmic behavior, such as early 
afterdepolarizations in a monolayer settings^\ Fourthly, optogenetic 
approaches offer new means to explore activity- dependent signaling 
in heart. Rapid electrical pacing can induce pathological changes 
leading to heart failure**^. Moreover, electrical resynchronization 
therapy can do much to revert myocytes from failing hearts towards 
a normalized profile characteristic of healthy hearts, according to a 
'reverse remodeling' process'*\ As well, activity- dependent signaling 
is believed to bias stem cells towards distinct cardiac lineages'*'^, and to 
coordinate development of embryonic and neonatal heart^^. The 
complex interaction of activity with multiple signaling pathways 
and genetic changes likely underlies these alterations in cellular 
phenotype^^"^". Optogenetic control of myocytes offers new flexibil- 
ity in specifying the activity profile of various preparations. Not only 
can excitable myocytes be altogether triggered or silenced, pre-excit- 
able cells can also be depolarized in controllable patterns, and the 
shape of action potentials in excitable cells adjusted. One intriguing 
signaling hypothesis is that changes in voltage itself (rather than 
ensuing alterations in calcium signaling) may trigger key pathways; 
optogenetic control during calcium channel blockade may prove 
revealing in this regard. Finally, optocurrents might ultimately be 
exploited as a next generation implantable antiarrhythmic device 
if means could be achieved to safely express optogenetic constructs in 
the heart. In this regard, it would be intriguing to drive expression of 
optogenetic constructs with promoters that are preferentially acti- 
vated in regions of pathology. Alternatively, injection of cardiac stem 
cells bearing optogenetic constructs could preferentially target 
regions of myocardial damage""^'^^. In all, the opportunities are vast 
and exciting, and should continue to draw creative minds to this 
field. 

Methods 

NRVM isolation. NRVMs were isolated as reported^^. Briefly, ventricles of P2 
Sprague-Dawley rats (Harlan) were excised and incubated (—16 hr, 4 C shaking 
water bath) in 50 mL HBSS (Gibco) containing 10 mg trypsin (-597,000 USP units/ 
g, Amersham Life Sciences). The next day, single ventricular myocytes were isolated 
by multiple rounds of mechanical shaking in 50 mL HBSS containing 50 mg 
collagenase Type 2 (215 units/mg, Worthington). Cells were filtered through a 70 
strainer (Falcon), and individual cells were re-suspended in M-199 medium (Gibco) 
containing 10% FBS (Gibco), 2 mM L-glutamine (Gibco), 19.4 mM glucose (Sigma), 
100 ]iM MEM non-essential amino acids (Gibco), 10 mM HEPES (Invitrogen), 
1.967 )iM Vitamin B12 (Sigma), and 12.84 [iM penicillin-G (1555 unit/mg, Sigma). 
The cells were successively plated in T75 and T175 flasks and incubated for 1 h each. 
Myocyte monolayers were cultured on plastic coverslips (McMaster, 85585kl03) pre- 
coated with 2.5 p-g/[iL fibronectin (Invitrogen), 2.5 X 10^ cells per 4.8-nmi diameter 
or 5 X lO'^ cells per 10-mm diameter coverslip. 

Molecular biology. For cardiac expression, eNpHR3.0 from pLenti-hSyn- 
eNpHR3.0-eYFP (from K. Deisseroth, Stanford University) was cloned into pPT- 
CMV-eGFP (from G. TomaseUi, Johns Hopkins University) and designated pPT- 
CMV-eNpHR3.0-eYFP. 

Lentivirus production and transduction. Lentiviruses were produced according to 
modified, published protocols'^^ using the calcium -phosphate method. For ChR2 
lentivirus production, HEK293T cells were transfected with 20 |j.g pLenti-EFlalpha- 
ChR2-YFP-WPRE (pLECYT), 10 ^ig deltaR8.74, and 10 ng VSVG (all from K. 
Deisseroth, Stanford University) per 15-cm plate. For eNpHR3.0 lentivirus 
production, 20 jig pPT-CMV-eNpHR3.0-eYFP, 7 |ig pMDLg-pRRE, 10 jig 
pMD2.G, and 5 |ag of pRSV-Rev (all from G. TomaseUi, Johns Hopkins University) 
were transfected per 15-cm plate. Virus was concentrated from supernatant using the 
ultracentrifuge method and stored in 15 uL aliquots at — 80"C. Titre was determined 
by transducing HEK293 cells in 10-fold serial dilutions and checking for fluorescence 
72 h later. For optical mapping, NRVM monolayers were transduced with either 
ChR2 (-4 i^L, 1 X 10^ pfu/mL) or eNpHR3.0 lentivirus (--5 |iL, 1 X lO^'^ pfu/mL). 
We did not find obvious signs that expression of either ChR2 or eNpHR3.0 affected 
resting potentials. In particular, action potential morphologies were largely 
unchanged (Figure 3D-F), arguing against appreciable depolarization resulting in 
decreased Na channel availability. Moreover, threshold voltages required for 
capturing preparations were unchanged with opsin expression (no opsin, 11.88 ± 
2.37 V (n — 30 coverslips); halorhodopsin, 11.74 ± 2.74 V (n = 35 coverslips); all 
n.s.), again suggesting maintained Na channel availability. 



Optical mapping. Optical mapping studies were performed on NRVM monolayers 
8-10 days post-plating, using di-4-ANEPPS^^ (Fig. 2b) or PGHl (all other mapping 
data). Coverslips of monolayers were placed in a custom-designed chamber and 
continuously superfused with warm (36 ± 0.5' C) Tyrode's solution (in mM): 135 
NaCl, 5.4 KCl, 1.8 CaCl2, 1 MgCl2-H20, 0.33 NaH2P04, 5 HEPES, 5 glucose, pH = 
7.4, 300 mOsm. Cells were stained with voltage -sensitive dye, and a Plexiglas cover 
placed over the chamber to minimize motion artifact from solution flow. Field 
stimulation was delivered by 2 parallel platinum line electrodes 2.7 cm apart. 
Transmembrane potentials visualized as optical fluorescent signals were collected 
directly below the chamber with an array of 253 optic fibers (1 mm diameter each) 
arranged in a 17-mm hexagon bundle, individually connected to photodiodes. 
Primary data acquisition and display were undertaken with software written in 
LabVIEW (Texas Instruments). GFP alone (or its color mutants fused to optogenetic 
constructs to visually track expression levels) were usually virally transduced in 
cultured NRVM preparations. As such, conduction velocities were generally slower 
than observed in untransduced preparations, likely due to a previously reported effect 
of GFP itself to diminish velocities'^'^. Conduction velocities may also have been 
slowed because we often enriched somewhat the usual fraction of cardiac fibroblasts 
in our NRVM cultures, to boost culture health and tolerance to strong expression of 
optogenetic constructs. 

For di-4-ANEPPS (Molecular Probes) mapping, cells were stained with 10 |iM of 
dye for 7 min. To record transmembrane potentials, the dye was excited perpen- 
dicular to the plane of the monolayer with an array of 26 high-power 530 nm LEDs 
(Kingbright) band-pass filtered at 530/50 nm. Fluorescence emission was long-pass 
filtered at 600 nm (PSCred; Brewer Science) before collection by the photodiode optic 
fiber array 

For PGHl mapping, cells were stained with 80 |J.M of dye (from Guy Salama) for 
9 min. To record transmembrane potentials, the dye was excited perpendicular to the 
plane of the monolayer with an array of seven high-power 655 nm LEDs (Luxeon 
Star) band-pass filtered at 690/60 nm (Omega Optical). Fluorescence emission was 
long-pass filtered at 760 nm (Newport, catalog no. 20CGA-760) before collection by 
the photodiode fiber optic array. For simultaneous excitation of PGHl and ChR2 or 
eNpHR3.0, the optical setup was modified for two orthogonal light paths using 
dichroic mirrors to split the red/blue or red/green lights at 475 nm (Omega Optical) 
or 618 nm (Omega Optical), respectively (Fig. 2a). ChR2 was excited with an array of 
seven high-power blue (447 nm) LEDs (Luxeon Star) band-pass filtered at 440/ 
60 nm (Omega Optical). eNpHR3.0 was excited with an array of seven high-power 
green (530 nm) LEDs (Luxeon Star) band-pass filtered at 560/50 nm (Omega 
Optical). The optics setup was controlled using a custom-built electronic circuit 
programmed for various light pulse protocols using MPLAB (Microchip MPLAB IDE 
version 8.76.00.00). 

When used. Bay K 8644 was diluted in Tyrodes to a final concentration of 1 |j,M Bay 
K 8644 and continuously superfused over the monolayers during transmembrane 
potential recordings. 

Whole-cell patch-damping. HEK293 cells were transduced with ChR2 or eNpHR3.0 
lentivirus and maintained as a stable cell line (non-monoclonal) for voltage-clamp 
recordings. Room temperature, whole-cell recordings were performed using an 
Axopatch 200A amplifier (Axon Instruments). P/8 voltage leak subtraction was used, 
with series resistances of 1-2 MQ and >70% compensation. Currents were lowpass 
filtered at 2 kHz (4-pole Bessel) and sampled at 10 kHz. The bath contained (in mM): 
138 NaCl, 4 KCl, 2 CaCU, 1 MgCl2-6H20, 0.33 NaH2P04, 10 HEPES, 5 glucose. The 
pipette solution contained (in mM): 129 K-gluconate, 10 HEPES, 10 KCl, 1 MgCl2, 0.5 
EGTA, 10 MgATP''. An array of seven blue (447 nm) LEDs (Luxeon Star) was used to 
excite ChR2 at a maximum power of 7.2 mW/mm^, and an array of seven green 
(530 nm) LEDs (Luxeon Star) was used to excite eNpHR3.0 at a maximum power of 
1.2 mW/mm^. A holding potential of 0 mV and an 8-s repetition interval were used. 
Data analysis via custom MATLAB software (Mathworks, MA). 

For cardiac APWs, optogenetic stimulatory light was turned on 500 ms before and 
after acquiring optogenetic currents during voltage-clamp with the APW. The APW 
utilized NRVM action potentials, as recorded using optical mapping with PGHl, and 
scaled as reported^^. 

Confocal imaging. NRVM monolayers bearing YFP-tagged ChR2 or eNpHR3.0 
were imaged with a 60X water objective on an Olympus FV300 confocal microscope. 
An Argon-515 nm laser line was used for excitation, and emission fluorescence was 
filtered at 535/30 nm. 

Data analysis of optical signals. Custom MATLAB software (Mathworks, MA) was 
used for post-processing of optical signals. AU signals were de-trended with a linear or 
polynomial fit and normalized to a range from 0 to 1. To select activation and 
repolarization times, optical signals were filtered using a low-pass Butterworth filter 
at 35 Hz. APD was calculated at 80% of repolarization (APDgo) (Fig. 5c). Average 
action potentials were generated by aligning waveforms at 50% of depolarization. 
Activation maps for di-4-ANEPPS and PGHl were constructed from activation times 
determined from the time of the stimulus to the instant at which waveforms first 
achieved 50% of peak. Colormaps represent activation times in ms. Activation maps 
for ChR2 triggered action potentials (Fig. 4a) were constructed from activation times, 
where time zero in the activation isochrone map corresponds to the start of the blue- 
light pulse. 
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Gap compensation for blue- and green-light bleed-through. At times, multiple 
pulses of blue or green light were delivered during the plateau and repolarization 
phase of the action potentials. The green- or blue-light pulses caused a baseline shift in 
our optical signals (Supplementary Fig. Si), as the photodiode array was positioned 
directly in line with illuminating light, and commercially available 760-nm, long-pass, 
emission filters imperfectly excluded the strong incoming light pulses. Accordingly, 
these periods of baseline shift were excluded and manifest as 'gaps' in our optical 
records. To estimate the trajectory of voltage during these gaps produced by blue-light 
pulses {Fig. 5), we utilized a custom curve-fitting algorithm. First, each action 
potential was linearly de-trended and normalized to a range from 0 to 1. Then, the 
control action potential (Fig. 5 (o)) was least- square-error fitted with a 
phenomenological function comprised of the weighted sum of three Hill equations: 



f"hi 



(1) 



Next, the elongation of APD due to a single blue pulse (Fig. 5c (i)) was fitted by a bi- 
exponential function from the time of the start of the blue pulse (Eq. 1). 



-Ffit.blLiepiUse ^ D- [exp( -{t- ipLLlse)/Tslow) - exp( -(t- ipulse)/Tfasl)] ^ 



(2) 



In particular, the difference of the action potential with and without the single blue 
pulse was least-square-error fitted with this bi- exponential function, and then this 
function (Eq. 2) was added to the control action-potential curve fit (Eq. 1). To obtain 
the curve fit of the two blue pulse action potential, the difference of the action 
potential with one and two blue pulses was fit with another bi- exponential function, 
and then this bi- exponential fit was added to the single blue pulse action potential 
curve fit (Fig. 5c (ii)). The same process was successively performed for three and four 
blue pulse action potentials (Fig. 5c (iii-iv)) until all fits were obtained (Fig. 5c, all 
curve fits on bottom right). In this manner, we minimized to a single gap the missing 
data to be estimated in fitting each pulse component (Eq. 2), thereby increasing 
precision. This 'gap compensation' procedure allowed for a continuous estimate of 
the voltage trajectory, despite the gaps apparent in optical voltage signals. 

For gaps induced by green-light pulses, the lower power of these pulses enabled a 
simple offset to be added to signals during baseline shifts, so as to produce a corrected 
signal that was approximately continuous at pulse onsets and offsets (Figs. 4 and 6). 
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